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ABSTRACT 

The chemical structure of neutral clouds in low metallicity environments is examined with 
particular emphasis on the H to H2 and C"^ to CO transitions. We observed near-IR H2 
(1,0) S(l), (2,1) S(l), and (5,3) 0(3) lines and the ^^CO J = 1 ^ line from 30 Doradus and 
N159/N160 in the Large Magellanic Cloud and from DEM S 16, DEM S 37, and LLSMC 36 in 
the Small Magellanic Cloud. We find that the H2 emission is UV-excited and that (weak) CO 
emission always exists (in our surveyed regions) toward positions where H2 and [C 11] emission 
have been detected. Using a PDR code and a radiative transfer code, we simulate the emission 
of line radiation from spherical clouds and from large planar clouds. Because the [C 11] emission 
and H2 emission arise on the surface of the cloud and the lines are optically thin, these lines are 
not affected by changes in the relative sizes of the neutral cloud and the CO bearing core, while 
the optically thick CO emission can be strongly affected. The sizes of clouds are estimated by 
measuring the deviation of CO emission strength from that predicted by a planar cloud model 
of a given size. The average cloud column density and therefore size increases as the metallicity 
decreases. Our result agrees with the photoionization regulated star formation theory by McKee 
(1989). 
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structure — Magellanic Clouds — radio lines: ISM 
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1. Introduction 

Stars form in dense, cold molecular clouds. Measuring the molecular gas content of the clouds is very 
important if we are to estimate the star formation efficiency and relate it to the properties of the clouds and 
to their environments. The total H2 mass, however, cannot be measured directly because the lowest levels 
of H2 from which the observable emission can arise have excitation energies (e.g., /S.E/k ~ 500 K, J = 
2^0) too high to be thermally excited in the cold (< 50 K) molecular clouds. In the Milky Way, the -'^^CO 
J = 1 ^ lincQ (hereafter CO J = 1 ^ 0) tr aces the molecular gas content. The conversion factor (X*^"^^) 
between the H2 column density and the velocity integrated intensity of CO has been measured via the virial 
theorem {X^"^^ = (2.5 - 8) x lO^o cm^V(Kkms-i), Solomon et al. 1987; Digel et al. 1997 and references 
therein), or via gamma-ray emission [X'^^'" = 1.1 — 2.8 x lO^*' cm^^/(Kkms^^), Bloemen et al. 1986; 
Digel et al. 1997 and references therein). The metallicity dependence of the conversion factor has been an 
issue. Cohen et al. (1988) and Wilson (1995) used cloud masses determined using the virial theorem to 
argue that the value of X increases as the metallicity of the individual galaxy decreases. Arimoto, Sofue, & 
Tsujimoto (1996) extend this conclusion to argue that there are radial increases in X in the Milky Way and 
M51 corresponding to radial decreases in metallicity. By contrast, Taylor, Kobulnicky, & Skillman (1996) 
showed that some low abundance galaxies have lower X, suggesting that a factor other than the abundance 
(e.g., temperature) can affect the measured value of X. 

Far-UV photons from massive young stars strike the surfaces of nearby molecular clouds^ and produce 
photodissociation regions or photon-dominated regions (hereafter PDRs, Tielens & HoUenbach 1985, 1997). 
In these surface layers, the far-UV photons dominate the ionization of atoms, the formation and destruction 
of molecules, and the heating of the gas. Inside the PDR, absorption by dust, C, and H2 diminishes the 
far-UV field. Several authors have constructed PDR models appropriate to conditions in the Magellanic 
clouds, with particular emphasis on the C+/C/CO transition (Maloney & Black 1988^; van Dishoeck & 
Black 1988b; Lequeux et al. 1994; Maloney & Wolfire 1997). In irregular galaxies, where metallicities and 
dust-to-gas ratios are lower than those in the Galaxy, far-UV photons penetrate deeper into clouds, and 
dissociate CO molecules to greater depths (Israel et al. 1986). Therefore, for a cloud with a given column 
density, the CO column density should be lower at lower metallicity. If the CO column density is high 
enough for the CO to self-shield against photodissociation {N{CO) > 10^^ cm^^, van Dishoeck & Black 
1988a), the CO column density will also be high enough for the CO J = 1 ^ line to be optically thick, 
and the CO J = 1 ^ line intensity (J Tmb dv) will not depend strongly on the metallicity. In that case, 
lower CO intensities can only stem from geometrical or beam-filling effects. On the other hand, if the cloud 
column density is not high, most of the CO will be dissociated and the resulting CO line will be optically 
thin and very weak. On the surface of the clouds, the destruction and formation of II2 molecules are also 
affected by the change of metallicity, but the mechanism is different from that for CO molecules. The II2 
molecules are dissociated by far-UV photons attenuated by dust or by II2 self-shielding. If II2 self-shielding 
dominates over dust attenuation, the H2 destruction rate is independent of the dust abundance. On the 
other hand, the H2 formation rate is proportional to the dust abundance, because H2 reforms on the 
surfaces of dust grains. 



The CO J = 1 — » intensity means the velocity integrated main beam brightness temperature, I = J Ta/s dv. 

^ In this paper, "a cloud" implies an UV-illuminated molecular structure either in isolation or in a complex. 

® Figure 4 in Maloney & Black (1988) is not correct. The figure should be replaced by Figure 2 in Israel (1988) or by Figure 
2 in van Dishoeck & Black (1988b). 
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The Magellanic Clouds are the best targets to test PDR models that include metallicity effects because 
of their proximity (d^^'~' = 50.1 kpc and d^^'^'^' = 60.3 kpc, Westerlund 1990), their low metal abundance 
(^ZLMC ^ 0.28, Z^*-^^ 0.54, = 0.050, and Z^^^^ = 0.21, where Z is normahzed to the Galactic 

value; Dufour 1984), and their low dust-to-gas ratio {p^^jl^ — 0.25 and p^.^l^ — 0.059, where p is normalized 
to the Galactic value; Koornneef 1984). In this paper, we observed the Magellanic Clouds in the near-IR 
H2 emission lines and in the CO J = 1 — > line (see Sections ^ and ||) . We compare the line intensities of 
H2 (1,0) S(l), CO J = 1 ^ 0, and [Cii] 158 ^m emission from the PDRs in the Magellanic Clouds with 
those from Galactic star formation regions (see Section U). Section ^ discusses the numerical PDR models 
which we compare to the observed data to learn how metallicity changes affect the chemical structure of 
the Galactic clouds and the clouds in the Magellanic Clouds. 



2. Observations 

Some limited regions in the Magellanic Clouds were previously observed in the H2 lines (Koornneef 
& Israel 1985; Israel & Koornneef 1988; Kawara, Nishida, & Taniguchi 1988; Israel & Koornneef 1992; 
Krabbe et al. 1991; Poglitsch et al. 1995). However, the published [Cii] and CO data (Johansson et at. 
1994; Poglitsch et al. 1995; Israel et al. 1996) cover more extended regions than the existing II2 maps. 
We observed near-IR H2 emission lines from the Magellanic Clouds with the University of Texas Near-IR 
Fabry-Perot Spectrometer whose equivalent-disk size^ {Oed = 81") is comparable to those of the existing 
[Cii] data {Oed = 68") and CO J = 1 ^ data {0ed = 54"). We also observed CO J = 1 ^ emission at 
positions where no emission had been detected at the sensitivity of the existing CO surveys. 



2.1. Near-IR H2 Emission Lines 

We observed the H2 (1,0) S(l) and (2,1) S(l) lines in 1994 December, and (1,0) S(l) and (5,3) 0(3) 
lines in 1995 October, at the Cerro Tololo Inter- American Observatory 1.5 m telescope using the University 
of Texas Near- Infrared Fabry-Perot Spectrometer (UT FPS, Luhman et al. 1995). The instrument was 
designed to observe extended, low surface brightness line emission, and has a single 1 mm diameter InSb 
detector to maximize the area-solid angle product. This product, Afi, depends on the telescope coupling 
optics but can be made as large as 4.5 x 10^'^ cm^ sr. The beam from the //30 secondary of the telescope is 
collimated and guided to an H-band or K-band Fabry-Perot interferometer and to a ~ 55 K (solid nitrogen 
temperature) dewar in which a doublet camera lens (focal length = 20 mm) focuses onto the detector. A 
0.5 % (for the H2 (5,3) 0(3) line) or a 1 % (for the H2 (1,0) S(l) and (2,1) S(l) lines) interference filter in 
the dewar selects a single order from the Fabry-Perot interferometer. We used a collimator of focal length 
— 838 mm in 1994 December run and a collimator of focal length — 686 mm in 1995 October. The change 
of the collimator affected the beam size and spectral resolution (Table ||) . Since the coupling of the 686 
mm collimator with the telescope optics is better, the beam profile in 1995 October is much closer to a 
box- function than the 1994 December beam profile (Figure 0). 

An automatic alignment routine aligned the Fabry-Perot etalon by executing every 5 — 8 minutes 
(Luhman et al. 1995). The Fabry-Perot etalon maintained alignment for 15 — 30 minutes, but the ambient 



^ The equivalent-disk size, Bed, is the diameter of a cylindrical beam whose solid angle is same as the integrated solid angle 
of the actual beam pattern. 
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temperature changes caused the plate separation to drift by the equivalent of 2 — 5 kms ^ per minute. 
Using telluric OH lines (Oliva and Origlia 1992), we calibrated the wavelength scale to within ±30 kms"^. 

We operated the Fabry-Perot interferometer in scanning mode at selected positions and in frequency 
switching mode at most positions. In the scanning mode, the plate separation of the etalon was varied to 
cover ±200 kms~^ centered at the H2 line in 15 sequential steps. Figure || shows observed H2 (1,0) S(l) 
and (5,3) 0(3) lines at the 30 Doradus (0, 0) position (see the object list in Table ||), and telluric OH 
lines. The OH (9,7) R2(2) (A = 2.12267 fim, wavelength in air, Oliva & Origha 1992) and the (9,7) Ri(l) 
(A = 2.12440 /im) lines are within the H2 (1,0) S(l) scan range. Figure || shows the red-wing of the OH 
(4,2) Pi(3) (A = 1.61242 ^m) line, and the OH (4,2) Pi(2) (A = 1.60264 ^m) line which was displaced by one 
free spectral range and penetrated through the blue side of the order sorting filter at Vlsr — 420 kms~^. 
The typical intensity of the OH (4,2) Pi (3) and (4,2) Pi (2) lines is 4 x 10"^ ergss"^ cm"^ sr~^ which is 
more than lO'^ times the H2 (5,3) 0(3) intensity. 

The OH intensity fluctuates spatially and temporally. The 1// power spectrum of the temporal 
fluctuations limits the sensitivity of the system, e.g., in K-band, the OH noise becomes important at 30 
seconds (coadded integration time) per step when the switching interval between the source position and 
an off-source "sky" position is two minutes (Luhman et al. 1995). In the observations of the Magellanic 
Clouds, we chopped the secondary mirror from the source to Aa = ±16' or —16' at 0.25 Hz in H-band and 
0.5 Hz in K-band. Each spectral step (at the same etalon plate separation) has one chopping cycle which 
consists of four exposures: object — > sky — > sky —^ object. Each exposure had an integration time of 0.5 or 
1 second. 

In the frequency switching mode, we tuned the Fabry-Perot interferometer to the line wavelength, Aon, 
and a nearby wavelength free of line emission, Xoff- One observing cycle consists of four steps: Aon — > ^off 

Xoff Aon- As is the case in the scanning mode, each step consists of one chopping cycle. We tried 
to place Xoff away from wings of the H2 instrument spectral profile for Aon and at positions away from 
telluric OH lines (see Table |). 

For flux calibration, we measured HR 1713 (B8I, mK = ±0.18 mag) in 1994 December and HR 8728 
(A3V, m^H = ±1.03 mag and mx — ±1.00 mag) in 1995 October. Even though the beam sizes of the 1994 
December run and the 1995 October run are different from each other, H2 (1,0) S(l) intensities measured 
at the same positions on both runs agree to within the errors. The absolute flux calibration is accurate to 
±20 %. 



2.2. CO J = 1 Emission Line 

We observed CO J = 1 ^ Une in 1995 December at the SESlf] located on La SiUa in Chile. The 
beam size (FWHM) of the SEST is 45". The CO intensities presented in this paper are the main beam 
brightness temperatures, Tmb (see the convention by Mangum 1993). The online temperature, T^, has 
been converted to T^/s by dividing by 0.7 (which is the product of the forward spillover and scattering 
efficiency, and the main beam efficiency). We used frequency switching to gain a factor of 2 in observing 
time and to avoid possible emission from reference positions. Since this method can leave residual ripples 
in the spectra, the detection limit for weak signals is determined by the ripple. For some positions, we 



The Swcdish-ESO Submillimeter Telescope, SEST, is operated jointly by ESO and the Swedish National Facility for Radio 
Astronomy, Onsala Space Observatory at Chalmers University of Technology. 
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complemented the frequency switched data with beam switched data with a reference beam about 12' away 
in azimuth. For example, in the 30 Doradus (0, —6') region, it turned out to be impossible to use the 
frequency switched data. 

We mapped CO J — 1 ^ spectra around positions where we detected H2 emission and where the CO 
emission was below the lowest contour level (3 Kkms~^) of previously published CO maps (Booth 1993; 
Johansson et al. 1994). The CO maps are fully sampled on 20" grids within each 81" H2 beam. In a single 
45" beam, the typical root-mean-square noise is 0.07 K for a channel spacing of 0.45 kms~^, and the typical 
la statistical uncertainty of the intensity I = J Tmb dv, integrated over a velocity range of 15 kms""'^, is 
0.2 Kkms"^ 



3. Results 

In low metallicity objects, the low dust abundance allows far-UV photons to penetrate for long 
distances beyond the Hii/H boundary. One possible scenario is that the transparency to UV photons leads 
to substantial regions of neutral gas where far-UV photons have dissociated CO but where self shielding 
permits hydrogen to be primarily in the form of II2. Especially in the N159/N160 complex, the published 
maps show the [C 11] 158 /xm emission (Israel et al. 1996) extending farther than the CO J = 1 — » emission 
(Johansson et al. 1994). These existing observations make such a possibility seem reasonable in the LMC. 

Near-IR H2 emission in response to far-UV radiation provides a direct test for the presence of molecular 
gas near cloud boundaries, albeit with no information about column density or total abundance. We 
selected observing targets in regions of the LMC where high spatial resolution V) [Cii] maps and CO 
J = 1 — * maps existed (Booth 1993; Johansson et al. 1994; Poglitsch et al. 1995; Israel et al. 1996). In 
the SMC, our II2 observations cover positions for which there were CO J = 1 ^ data available (Rubio et 
al. 1993) but do not coincide with the observed [Cii] positions (Israel & Maloney 1993). Table || lists the 
observed sources and their reference (0, 0) positions. 



3.1. H2, CO, [Cii], and far-IR 

The II2 (1,0) S(l) line intensity in the surveyed regions is < 4 x 10~^ ergss"^ cm~^ sr~^, except in the 
central regions of 30 Doradus and N160. About 70% of the observed points have detections of the (1,0) S(l) 
line with a significance of 2a or more. The II2 (1,0) S(l) intensities are listed in Table ^, and the II2 
(2,1) S(l) and (5,3) 0(3) intensities in Table ^. Reddening toward the stars in the LMC has been measured 
using (U — B) and {B — V) colors. The Galactic foreground extinction. Ay, toward the LMC is 0.23 ± 0.07 
mag (Greve, van Genderen, & Laval 1990; Lee 1991). Applying Ak = 0.112 Av(Rieke & Lebofsky 1985), 
the foreground Ay implies a negligible Ak of ~ 0.03 mag. The molecular clouds in the LMC may, however. 



extinguish the H2 emission from their own back sides. Section 5.3 discusses the II2 emission from the back 
side of the clouds. 

Table I hsts the CO J = 1 ^ and [C 11] intensities convolved to the beam size of the II2 observations 
{Sed ~ 81"). We assumed that the beam profile of the University of Texas Fabry-Perot Spectrometer is a 
box function. 

The columns for far-IR intensity and dust temperature in Table are from the IRAS data processed 
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at IPACf] using standard HIRES processing. The HIRES processor enhances the spatial resolution of IRAS 
images using a Maximum Correlation Method algorithm (Aumann, Fowler, & Melnyk 1990). The effective 
resolutions are l.'O — 1'.6 at 60 /im and 1'.6 — 2.'4 at 100 /zm. We calculate the integrated far-IR intensity 
using the approximation of Helou et al. (1988): 

Ifir = 1.26 X 10"^ (2.58 + /lOO^m) : (1) 

where /eo/^m and /looMm are in units of MJy/sr and Ifir is in ergss"^ cm~^ sr~^. This equation is valid 
within 20 % errors, when 31 < Tdust < 58 K. We deduce the dust temperature, Tdust, from the ratio of 60 
fim to 100 /im intensity: 

leafim _ f J^60/jm_ \ exp(feiyino^m/fcJdnst) ~ 1 

where h is the Planck's constant, k is the Boltzmann's constant, and v is the frequency in units of Hz. 
In the above equation, we assumed that the 60 /zm emission is not dominated by small (D c± 5 nm), 
thermally-spiking, grains but by large grains in steady state temperatures (Draine & Lee 1984), and that 
the dust emissivity (Q^) is proportional to i/" where the dust emissivity index, n, is 1. 



3.2. H2 Results 

3.2.1. SMC 

In the SMC, we detected H2 (1,0) S(l) emission, with S/N > 2, near IRAS sources in the DEM S 
16, DEM S 37, and LI-SMC 36 regions. Toward supernova remnants, e.g., DEM S 37 (+0.'38, -3f35) and 
LI-SMC 36 (—3', 0), we did not detect H2 emission at a la level of 2 x 10~^ ergss^^ cm~^ sr~^. 



3.2.2. LMC: 30 Doradus 

Figure H shows the far-IR, CO J = 1 ^ 0, and H2 (1,0) S(l) intensity maps of the 30 Doradus region. 
At the 30 Doradus (0,0) point, the intensity of the H2 (1,0) S(l) line, Ih^, is 1-1 x 10~^ ergss"^ cm~^ sr~^ 
in our 81" beam (hereafter, denotes the intensity of the H2 (1,0) S(l) line). Poglitsch et al. (1995), 
using their imaging NIR-spectrometer FAST, observed intense H2 (1,0) S(l) emission around the central 
cluster R136 and showed that the H2 source appears highly fragmented (< 5" or 1 pc scale) with a typical 
intensity of ^ 1.6 x 10~^ ergss"^ cm~^ sr~^. 

Our observations show that the H2 emission in 30 Doradus is very extended, > 40 pc, compared to the 
known extent of the CO J=1^0 emission (~ 35 pc, see the right map in Figure H). At (or 22 pc) from 
the (0,0) position, the H2 (1,0) S(l) intensity is only a factor of two lower than at the peak (see Figure ||). 
We also detected faint H2 (1,0) S(l) emission at (0,-6'), = 1-6 x 10"^ ergss^-^ cm~^ sr~^, where CO 
J = 1 ^ emission had not been detected during the survey of the ESO-SEST Key Program (Booth 1993). 



^ IPAC is funded by NASA as part of the IRAS extended mission program under contract to JPL. 
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3.2.3. LMC: Nl 59/1 60 

The N159/N160 Hii complexes are 40' south of 30 Doradus. The [Cii] hne (Israel et al. 1996) and 
the far-IR continuum distributions (left map in Figure |^) show that the far-UV fields are strong both near 
N159 and N160. On the other hand, the CO intensity around N159 is more than four times stronger than 
that around N160. 

The H2 (1,0) S(l) line had been detected at "N159 Blob" (a compact Hii source with a size of 8" x 6", 
Heydari-Malayeri & Testor 1985) by several groups. In Table ^, we compare the previous data with our 
new results. The flux increases as the beam size increases. This suggests that, if there is a single source, 
the emission region is more extended than ~ 20". Alternatively, there could be clumpy emission filling our 
81" beam with an area covering factor ~ 20% of the covering factor in the inner 6" x 6" region. In our II2 
survey, we observed very extended (> 5' or 70 pc) II2 (1,0) S(l) emission from the N159/N160 Hii complex 
(see the right map in Figure |). In the N159 region, we detected H2 emission where the CO cloud complex 
is bright and extended (see Figure 1 in Johansson et al. 1994). In the N160 region, however, we also 
detected II2 emission beyond the lowest CO J = 1 ^ contour level (3 Kkms~^) in the map by Johansson 
et al. (1994)1^. In spite of the weak or absent CO emission in the N160 region, the II2 observations indicate 
that the size of the molecular cloud complex is as big as that in the N159 region. 

3.3. CO J = 1 ^ 

At several positions in the outer regions of 30 Doradus and N160, we detected II2 emission where 
earlier CO surveys failed to detect the J = 1 ^ line. In order to determine if all CO is dissociated at these 
positions, we observed the regions again in the CO J = 1 ^ line (see Section |2^ ) with higher sensitivity 
{a ~ 0.2 Kkms"^) than the previous observations of which the lowest contour level was 3 Kkms~^ (Booth 
1993; Johansson et al. 1994). 

We made a fully sampled CO map inside the UT FPS beam at (0, —6') in 30 Doradus, and detected 
CO J = 1 ^ at a level of 0.5 Kkms"^ (see Figure ||). We also made CO maps in the outer regions of 
N160 where the previous CO J = 1 ^ map did not show any CO emission, and sampled some positions 
in N159 to confirm the flux calibration of the new observations (see Figures ^ and||). We plot the contours 
at logarithmic intervals to emphasize the edges of the molecular cloud complexes; note the dense contour 
lines on the northwestern side and on the western side of the N160 complex. The CO J — 1 emission 
regions cover the H2 emission regions except the (— 2f5, +7') and the (—2', +10') positions where the CO 
cloud complex fills 50 — 70 % of the II2 beams. 

3.4. H2 Excitation Mechanism 

Table I lists the observed H2 (2,1) S(l)/(1,0) S(l) and (5,3) O(3)/(l,0) S(l) line ratios in the 
Magellanic Clouds. Observations of the H2 lines in archetypal shocked regions, e.g., Orion BN-KL, HH 7, 
and the supernova remnant IC 443, show that the H2 ratios of (2,1) S(l)/(1,0) S(l) are almost constant 
~ 0.08 or Texc ~ 2 X 10^ K (Burton et al. 1989; Richter, Graham, & Wright 1995). Assuming that the 



See also Figure 2 in Israel et al. (1996) which used the data of Johansson et al. (1994). More recently, the SEST Key 
Programme has observed more positions in 30 Doradus and N159/N160 (Johansson et al. 1998). 
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excited levels are in LTE and T^^c = 2000 K, the (5,3) O(3)/(l,0) S(l) ratio should be only - 9 x 10"^ 

Molecular hydrogen in PDRs absorbs 91.2 — 110.8 nni photons in the -B^Sj — X^'^'^ Lyman and 
C^n„ — ^^^g Werner bands. About 15% of the electronically excited molecules are dissociated (Draine & 
Bertoldi 1996). The remaining 85% of the excitations result in populations of various ro- vibration levels 
of the ground electronic state. If n(H2) < 5 x 10^ cm^"^, the relative line intensities arising in UV-excited 
H2 are insensitive to density or to UV field strength (Black & van Dishoeck 1987). In this pure- fluorescent 
transition case, the H2 ratio of (2,1) S(l)/(1,0) S(l) is 0.56 and (5,3) O(3)/(l,0) S(l) is 0.38. At densities 
> 5 X 10^ cm~^ (Luhman et al. 1997), the coUisional de-excitation of UV-pumped II2 begins to affect the 
ro-vibrational level populations. If the PDR boundaries become sufficiently warm [T > 1000 K) and dense 
(> 5 X 10* cm^"^), coUisional excitation thermalizes the low vibrational {v < 2) level populations (Sternberg 
& Dalgarno 1989; Luhman et al. 1997). 

The detections of II2 (5,3) 0(3) line in the Magellanic Clouds (Table ^) verify the UV-excitation of 
II2. The observed H2 line ratios from the N160 region are those expected for pure- fluorescence. The ratios 
toward other regions show that the II2 ro-vibration levels may be affected somewhat by collisions. In 30 
Doradus, the peak H2 (1,0) S(l) intensity is 2.3 x 10"* ergss"^ cm~^ sr~^ (Poglitsch et al. 1995) which 
is brighter than the maximum predicted by our PDR models (1^2 ^ 9.6 x 10~^ ergss"^ cm~^ sr~^, see 
Figure [l^) . The models may underestimate the intensity by as much as a factor of two in regions with high 
densities and high UV fields by neglecting the large increase in UV pumping in the very warm (420-600 
K) region at the cloud edge. Alternatively, there may be some coUisional excitation of H2 v=l as well 
as coUisional de-excitation of the high v states, as one finds in clouds subjected to intense far-UV fields 
(Luhman et al. 1997). 



4. Comparing v^rith Galactic Clouds 

We compare the H2 (1,0) S(l), [Cii] 158 /im, CO J = 1 ^ 0, and far-IR emission from the clouds in 
the Magellanic Clouds (see Table ^ to the emission from star forming clouds in Orion and NGC 2024 for 
which we have comparable data sets (see Table 0). In the LMC, we select positions where complete II2, 
[Cii], and CO data sets exist: 5 positions in the 30 Doradus region, 4 positions in the N159 region, and 8 
positions in the N160 region. In the SMC, we use 4 positions with only CO and II2 data. 



4.1. Data from NGC 2024 and Orion A Star Formation Regions 

For the Galactic cloud data, we make use of the published II2 (1,0) S(l), [Cii] 158 /im and CO 
J 1 data in Schloerb & Loren (1982), Stacey et al. (1993), Jaffe et al. (1994), Luhman et al. (1994), 
Luhman & Jaffe (1996), and Luhman et al. (1997). The far-IR data are from the HIRES processed IRAS 
data (see Section ^.l[ ). In the Orion molecular cloud, Stacey et al. (1993) made two [Cii] strip maps in 
Right Ascension, both of which pass close to the CO J = 1 ^ peak. The [C 11] flux along the strips 
was "bootstrapped" , integrated by assuming zero flux at the ends of the cut and summing the chopped 
differences. The data from cut 1 (observed west to east) and cut 2 (observed east to west) are not in 
complete agreement; therefore, we take only those data with: 

^UcvAl^Iau^ < 0.4 . (3) 
Icutl + Icut2 
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We also exclude positions where the IRAS data are saturated, i.e., leo^m > 1500 MJy/sr, and where the 
H2 ro-vibrational level populations begin to show effects of coUisional de-excitation ((6,4) Q(l)/(1,0) S(l) 
< 0.15, see Table 2 in Luhman & Jaffe 1996). Some positions where coUisional de-excitation was unlikely 
were not observed in the H2 (6,4) Q(l) line but are included in the data we use (see the discussion in 
Section 5.5.1). Table |^ lists the compiled data sets in the Galactic Clouds: 4 positions from the Orion A 



molecular cloud (hereafter "Orion"), and 11 positions from the cloud associated with the NGC 2024 Hii 
region (hereafter "NGC 2024"). 

Figure ^ compares the data from 15 positions in the Galactic clouds with the data by Jaffe et al. 
(1994) who presented CO J = 1 ^ and [C 11] data from NGC 2024, and showed that the different parts 
of the source have very different distributions on a Icii versus Ico plot. The left plot in Figure ^ shows 
the distinction between the cloud proper zone {Aa > —10' with respect to NGC 2024 IRS 1, shown as open 
circles) and the western edge zone (Aa < —10', shown as open triangles) in NGC 2024. Data from the 
cloud proper zone agree with the PDR models of Wolfire et al. (1989), while many of the Icii/Ico ratios 
toward the western edge zone are much higher than any of the model ratios. Jaffe et al. (1994) interpreted 
the western edge zone results as implying that the mean column density of clouds decreases to the west. 
The right plot in Figure || shows the location in the Icii — Ico space of the data sets in Table |^, which will 
be used to compare with those from the Magellanic Clouds. 



4.2. Far-IR Data 

4-.2.1. luv versus I fib. 

The intensity of the interstellar far-UV radiation field is usually expressed in terms of the scaling 
factor luv , the mean intensity in the solar neighborhood (Draine 1978; also see footnote 2 in Black & van 
Dishoeck 1987). The critical part of the UV range for H2 fluorescence, C ionization, and CO dissociation 
is 91.2 < A < 113 nm (Black & van Dishoeck 1987; van Dishoeck & Black 1988a). The intensity of the 
luv = 1 field integrated over 91.2 < A < 113 nm is 3.71 x 10^^ ergss~^ cm~^ sr~^. Note that some authors 
(e.g., Tielens & HoUenbach 1985) use Habing's (1968) determination of the local far-UV field, 1.3 x 10"'' 
ergss"^ cm~^ sr~^, for the integrated intensity between Q < hv < 13.6 eV or 91.2 < A < 206.6 nm, and 
use the symbol Go to indicate the degree of enhancement over this standard integrated intensity. The 
corresponding intensity of the Draine field integrated over the same band as used for Go is 2.13 x 10"'* 
ergss"^ cm~^ sr~* (for a detailed comparison between the Draine field and the Habing field, see Draine & 
Bertoldi 1996). 

Most of the far-UV energy is absorbed by grains and reradiated in the far-IR: 

Ifib. = 2 X (2.13 x 10"'') luv ergss'^ cm^^ sr"!, (4) 

where the factor of 2 accounts for incident radiation at longer wavelengths, /^(A > 206.6 nm) (Wolfire, 
HoUenbach, & Tielens 1989). 



4-2.2. Ifir versus Tdust 

We can use the far-IR emission to normalize the various line intensities to compensate for beam filling 
factor effects. Before we do so, however, we need to understand how the far-IR emission arises. Figure 
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shows Ipi versus Tdust for the Galactic star forming cloud positions and the Magellanic Clouds positions in 
our sample. We derive Ifir f^rom the measured 60 /im and 100 /im intensities using Equation ^ and derive 
Tdust using Equation ||. Assuming that the dust is heated by an external far-UV field, we can approximate 
the results of PDR models: 

Tdust = 13.5 {luv 3§^) K (5) 

where Te// is the equivalent stellar surface temperature which would produce the incident UV field 
(HoUenbach, Takahashi, & Tielens 1991; Spaans et al. 1994). Assuming Tf.ff = 3 x 10"* K, the expected 
relation between Tdust and Ifir is: 

log/F/fl - -9.02 + h\ogTdust, (6) 
which is shown as a dotted line in Figure 0. 

The observed Tdust and Ifir distributions in Orion and NGC 2024 agree with the model in Equation ^. 
Ifir is proportional to the beam filling factor, while Tdust, which was deduced from the ratio of 
I m lira 1 1 wo ^va, is independent of the beam filling factor. With a beam size of ^ 1' (or ~ 0.12 pc at the 
distance of Orion; 415 pc, Anthony-Twarog 1982), the projected beam filling factor for the clouds in Orion 
and NGC 2024 is ^ 1, which explains the agreement between the observed Tdust — Ifir relation and the 
model. If we assume that the dust size distribution is independent of metallicity and that the clouds are 
optically thick in the far-UV, the Tdust-lFiR relation should be independent of the dust abundance (see 



Section |5.2| for more discussion). For the Magellanic Clouds, the observed Ifir values are an order of 
magnitude weaker than the values predicted by the Tdust — Ifir relation. This difference implies a beam 
fining factor in the Magellanic Clouds of - 0.1 (d^*^"^ = 50.1 kpc, Westerlund 1990). 



4.3. Line Intensities Divided by Ifir 

The line and continuum emission we describe here arises in the layers of molecular clouds where 
UV photons influence the chemistry and the physical conditions. Various parameters affect the emergent 
intensity: 

li = -qei fi{nH,Iuv,SvD,Z), (7) 
Ifir = t] eFiRi.26 x 10^^ luv , 

where li is the observed CO J = 1 — > 0, [Cii] 158 /im, or H2 (1,0) S(l) emission line intensity, and Ifir. is 
the observed far-IR continuum intensity in units of ergss"^ cm~^ sr^^. fi is the combined intensity arising 
from the front and back sides of a plane-parallel model cloud which fills the beam. The front and back 
surfaces of the cloud are exposed to external far-UV radiation, and are perpendicular to the line-of-sight. fi 
depends on the hydrogen number density, hh — ri(H) -I- 2n(H2), where ri(H) and ni^-i) are the atomic and 
molecular hydrogen density, respectively; the far-UV field strength, Iuv\ the Doppler velocity dispersion, 
8vf) (or half- width at 1/e point); and the metal abundance, Z. r\ is a beam filling factor and e is a geometric 
correction factor which we describe below. 

Unlike the substructure of the fully molecular interior of clouds, which may have many size scales, the 
penetration scale length of far-UV photons and UV photochemistry insure that structures within giant 
molecular cloud (CMC) complexes with UV-illuminated surfaces and CO-bearing cores have a minimum size 
determined by their density and metallicity. Measurements and theoretical arguments show PDR structures 
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outside the dense cores of molecular clouds are clumpy on a size scale of ~ 1 pc (Burton, HoUenbach, & 
Tielens 1990; Jaffe et al. 1994), comparable to or larger than the Orion beam but very much smaller than 
the beams in the LMC. Since the measured intensity is a beam average over any source structure, we have 
to correct for the effects of different beam filling factors in order to directly compare the emitted intensities 
from the LMC with those from the Galactic clouds. We define the beam filling factor, ry, as the fraction of 
the observed beam area filled by a single cloud (not an ensemble of clouds), using the outermost edge as the 
cloud boundary (or the boundary between Hii and Hi regions): 

(2Rcioud\^ 

where Rdoud is the radius of the cloud to the outer (Hii/H) boundary, d is the distance to the cloud, and 
Bed is the diameter of the telescope equivalent disk in radians, r] is same for all types of emission in the 
same cloud. In order to simplify model simulations in Section |^, we assume that only one cloud is in the 
telescope beam, the telescope has a box-function beam profile, and the telescope beam area is always larger 
than the cloud size, i.e., d Bed ^ "^Rcioud or 77 < 1. 

It may be more realistic to consider three dimensional rather than planar clouds. In case of a three 
dimensional cloud with an external far-UV field whose intensity is uniform on the surface of the cloud, 
the far-IR, [Cii], and II2 emission arise in the outer shells of the cloud, while the CO J = 1 ^ emission 
arises from the surface of the CO core inside the cloud. We use a "geometric correction parameter" , or 
tpiR, in Equation ^ to simulate the observed intensity of a three dimensional cloud. The spherical geometry 
parameter accounts for limb-brightening and for the size differences between the CO, [Cii], H2, and far-IR 
emission regions in a given spherical cloud. 

When examining our three dimensional model clouds, we can eliminate the beam filling factor, ry, by 
dividing li by Ifib.'- 

Ij _ hjnH, Iuv,Svd,Z) 

Ifir ~ epiR 4.26 x IQ- ^ luv ' 
We will therefore divide the observed intensities by the far-IR intensity at each position as a means of 
removing distance-related beam filling factor effects in our subsequent analysis of the data from the Galaxy 
and the Magellanic Clouds. 



4.4. Relationship between /co? Iciit and Ih2 

Table ||shows the average ratios between the observed H2 (1,0) S(l), CO J = 1 ^ 0, and Cii intensities, 
and the standard deviations of the ratio distributions in the Galactic clouds and in the Magellanic Clouds. 

Figures |, |, and |l^ show plots of IH2/IFIB. versus Icii/Ifir-, Icii/Ifir versus Ico/Ifir, and 
Ih-i/Ifir versus Ico/Ifir, for the Galactic clouds and for the clouds in the Magellanic Clouds. From these 
figures and Table ||, the line ratios of \og[lH2/ Icii) between the Galaxy and the Magellanic Clouds are in 
good agreement, and the ratios of log(/c////co) and log(///2 //co) in the Magellanic Clouds are shghtly 
higher (but in agreement within the standard deviations) than those in the Galaxy. Section |^ will discuss 
our PDR models and compare the observed data with PDR models with various input parameters. 
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5. Models 

Using plane-parallel codes, the metallicity dependence of PDR structure and of emergent line intensities 
has been calculated and discussed by several authors (Maloney & Black 1988; van Dishoeck & Black 
1988b; Wolfire, HoUenbach, & Tielens 1989; Maloney & Wolfire 1997). The emergent intensities of the 
H2 (1,0) S(l), [Cii], and CO J = 1 ^ lines from PDRs depend on nn, luv, Svd, and Z (see Section 



4.3). As long as the CO J = 1 ^ line is optically thick, the resulting intensities of CO, [Cii], and H2 
(1,0) S(l) lines are not very sensitive to the metallicity. A simplified analysis can illustrate the reason 
for this insensitivity. In the outer part of the PDR, gas-phase carbon is in the form of C"*". When the 
metallicity is lowered, the number density of C"*" ions drops, but the far-UV photons penetrate deeper 
into the clouds because the dust-to-gas ratio is also lower. Thus, the column density of C+ in the PDR 
is almost independent of the metallicity. The [C 11] line is optically thin, and the intensity is proportional 
to the column density of 0+. The CO J = 1 ^ intensity depends much more weakly on the CO column 
density because the line is optically thick and most of emission arises on the surface of the CO region inside 
the cloud. Therefore, the CO J = 1 ^ intensity outside of very high column density cloud cores does not 
depend strongly on the metallicity when one considers emission from a plane-parallel slab. We will discuss 



the metallicity dependence of H2 in Sections 5.2 and 5.3 



As we discussed in Section O, in the case of three dimensional clouds, we use the parameter, or 
^FiR, to account for effects such as limb-brightening and the size differences between the outer shells from 
which far-IR, [Cii], and H2 (1,0) S(l) emission arises and the inner CO core from which CO J = 1 ^ 
emission arises. The depths from the cloud surface to the C+/C transition layer and to the H/H2 transition 
layer are about inversely proportional to the dust abundance, so the metal and dust abundance are more 
important parameters in spherical-shell PDRs than in plane-parallel models. Storzer, Stutzki, & Sternberg 
(1996) and Mochizuki (1996) modeled a PDR on the surface of a spherical cloud with Galactic metallicity. 
This Section presents our own models of spherical-shell PDRs and calculate the emergent line intensities for 
a range of densities, UV fields, and metallicities. 



5.1. Codes 

We first ran a plane-parallel PDR code (van Dishoeck & Black 1986; Black & van Dishoeck 1987; 
van Dishoeck & Black 1988a; Jansen et al. 1995) with a range of densities {nn — 5 x 10^, 5 x 10'^, and 
5 X 10* cm-3), UV fields {luv = 10, 10^ 10^, and 10"*), and metallicities (for the Galaxy, the LMC, and 
the SMC). In this code, one side of the model cloud is exposed to UV radiation, and the cloud is divided 
into 200 slabs, each of which is in chemical steady state. Since the code only includes inelastic collisions 
of II2 within v=0, the ro-vibrational level populations of H2 are not correct at uh > 5 x 10'' cm~'^. The 
PDR code calculates the gas temperature and chemical abundances in each slab. Since there is no detailed 
information on the grain properties in the LMC and SMC, we assume that the formation rate of H2 on 



grains scales linearly with the dust abundance (section 5.2) and that the heating efficiencies of the grains 
are the same as in the Galactic case. In practice, the lower H2 grain formation rate is accomplished in the 
models by decreasing the parameter yf (Black & van Dishoeck 1987) from 1.0 to 0.25 and 0.1 in the LMC 
and SMC, respectively. Table ^ lists input parameters for the code. 

We employ spherically symmetric cloud models to study the effects arising in three dimensional chemical 
structures. The outputs, e.g., chemical abundances and kinetic temperatures, from the plane-parallel 
code are applied to the spherical shells of the cloud used in the radiative transfer model. We map the 
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temperatures and abundances derived at each distance from the Hii/Hi interface by the plane-parallel 
model into the spherical radiation transfer model, ignoring any changes in the chemistry or thermal balance 
due to the difference in geometry. The Monte Carlo code of Choi et al. (1995; hereafter MC) calculates 
the level populations of the atoms and molecules. The MC code simulates photons in a one- dimensional 
(spherical) cloud, and adjusts the level populations according to the result of simulations until the 
populations converge. We assume a purely turbulent velocity field with a Doppler velocity dispersion of 1 
kms~^ {6vd or half- width at 1/e point) with no systematic motion. Since the MC calculation only includes 
40 slabs, the 200 slab plane parallel model was smoothed in such a way as to retain high resolution at the 
transition regions (Li 1997). 

Wc use the output of MC to calculate emission lino profiles using a Virtual Telescope code (Choi et 
al. 1995; hereafter VT). The VT code convolves the integrated emission from each spherical shell along 
the line-of-sight with a virtual telescope beam profile to simulate observations. While it would be better to 
examine the geometric effects using fully self-consistent models which calculate the detailed H2 excitation, 

CO photodissociation, and CO and Cii line intensities, our present three-step procedure should be good 
enough to allow us to analyze the global properties of Galactic and LMC clouds and to establish relative 
trends. 



5.2. Effect of H2 Self-Shielding 

Inside neutral clouds, the external far-UV field is attenuated by dust absorption, and by C and H2 
absorption. CO absorption of far-UV is negligible in the outer part of the cloud because most of the 
carbon is ionized. H2 can survive in the outer parts of the cloud, either as a result of self shielding or as 
a result of shielding by dust. We can analyze the conditions under which H2 self-shielding from far-UV 

photons is dominant over shielding by dust (adapted from Burton, HoUenbach. & Ticlens 1990). Within a 
plane-parallel cloud, the H2 formation rate, F{x), at depth x measured from the surface of the cloud toward 
the center is 

F{x) =n{R,x)nHq^ , (10) 

where uh is the hydrogen number density, nn = n(H, x) -|- 2n{B.2,x), assumed to be constant over the 
cloud; n(H, x) and n(H2, x) are number densities, at x, of H atoms and H2 molecules respectively; K is the 
value of the H2 formation rate coefficient (a quantity which depends linearly on the dust abundance) for 
the Galactic dust abundance; q is the dust abundance normalized to the Galactic value, i.e., g'^^^ = 1. 

is a slowly varying function of the gas temperature (SR oc T^/^), so we take an average value (3 x 10"-'^'^ 
cm^'^ s^^. Burton, HoUenbach, & Tielens 1990) for this analysis. If we assume that, at the position we 
are considering, the dust optical depth is negligible and the H2 absorption is governed by the square root 
portion of the curve of growth (Jura 1974) the far-UV field is attenuated by A''(H2, a;)^/^, where 



N{R2,x) = / n{ll2,x)dx . (11) 
Jo 

The H2 destruction rate, D{x), at depth x is 

where lo is the unshielded dissociation rate of H2 at luv = 1 (7.5 x 10~^^ s~^. Black & van Dishoeck 
1987); and /? is the self-shielding parameter (4.2 x 10^ cm~^, Jura 1974). We can integrate the steady state 
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equation, F(x) = D{x), over x: 

nnq^ [ n(H, x) dx ^ luv lo (3 f "(^2, a:) ^ 

By substituting Equation |ll| into the above equation, we get: 

nHq^N{R,x) ^2Iuv Iof3N{R2,xy/^ . (14) 

H2 self-shielding is more important than dust shielding, when the far-UV field attenuation by dust is still 
negligible, 

Tdust{Xo) < - , (15) 

at the point, Xo, where the molecular hydrogen column density becomes equal to the atomic hydrogen 
column density, 

iV(H,Xo) = 27V(H2,Xo) . (16) 

Tdust{x) is the optical depth of dust at A = 100 nm. Even though the ratio of Tdustf-^v depends on the 
chemical composition and the size distribution of the dust, the change of this ratio from the Galaxy to the 
LMC and the SMC is negligible compared to that of the dust-to-gas ratio, e.g., the ratios of Ao.i^m/AB in 
the Galaxy, the LMC, and the SMC are 3.4, 4.1, and 5.1, respectively (Pei 1992), while Ay/Nn changes by 
a factor of 10 (Table 9). To simplify our analysis, we assume that the ratio of Tdust/Ay is constant and that 
only the dust-to-gas ratio (adopted from Bohlin et al. 1983 and Black & van Dishoeck 1987) depends on 
the metallicity: 

Trf„,t(x) ^3.0Av(x) , (17) 

and 

Av{x) = 6.29 X 10-^^pdustNH{x) , (18) 

where pdust is the dust-to-gas ratio normalized to the Galactic value (see Table 9), and Nh{x) is the 
hydrogen column density, Nh{x) ~ N(H,x) + 2N(ll2,x). Substituting Equation ^ into Equation Q we 
obtain the following relationship between the column density at which self shielding becomes effective, and 
the density, metallicity, and strength of the incident UV field: 

M qnn 

~ 2.1x10^2:^. (19) 
quH 

From Equations |l^, |l^, |8|, and |l^, we derive the conditions under which the H2 self-shielding is dominant 
over shielding by dust: 

^ > 1.3 X lOV.-i^/ , (20) 

where we assume that the optical absorption properties and the efficiency for H2 formation of the dust per 
unit H atom vary in the same way as with dust abundance: q = pdust- Figure |l^ shows the linear size of 
the C"*" region {X(j+) and the region {X^*: H2 denotes vibrationally excited H2) from the FDR model 
results: Xc+ and Xfj* are the distances from the surface of the cloud to the inner edges of the C"*" region 
where n{G'^ , Xq+) — n{CO,X(j+), and of the H2 region where 7i(H, AT/fj) = 7i(H2,Xh*) respectively. 

The plot of Xh* in Figure |l^ shows how H2 self-shielding affects the depth of the region from the 
surface of the cloud. Based on Equation EG, the dotted lines divide the {Ijjv, nn, Xh' ) space into a region 
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where dust absorption is dominant and a H2 self-shielding dominant region. We can predict the behavior of 
Xh^ as a function of luv, nn, and pdust, using the relations derived in this section. When dust absorption 
is more important than H2 self-shielding, the H2 zone ends when the far-UV field is attenuated to a fixed 
level, independent of the incident field. As a result, changes in Ijjy result in a variation of r^ust according 
to: 

luv exp(-Trf„st) = C , (21) 
where C is a constant. Substituting Equations [l^ and |l^ into the above equation gives that 

NniXn;) ^ 1.22 x lO^i log^^vlogg (22) 

Pdust 

where Nh{Xh*) is the hydrogen column density at Xh*- In other words, if we increase luv by an order of 
magnitude, Nh{Xh') increases only by an additional additive factor of 1.22 x 10^^ Pdust cm"^. 

When H2 self-shielding is dominant, from Equation pj|, 

logiVH(Xffj)~24.6 + 21ogf^') , (23) 

\qnH J 

e.g., if luv decreases by an order of magnitude, Nh{Xh*) decreases by two orders of magnitude. The 
above equation explains, why, when H2 self-shielding dominates, the depth of the H2 shell decreases rapidly 
as luv decreases or as uh increases (see Figure pd] ). 

In the LMC, (Pduif )"^^^ - 2-2 (see Equation |2^), and the dotted line is shifted toward lower luv 
by only a factor of 2.2. Even in the case of the SMC, {Pdusf)~^^^ — 3-2. Variations in metal abundance, 
therefore, do not significantly affect the H2 self-shielding criterion. 



5.3. Emission Intensities without Spherical Geometry Effects, e 

We first ran the MC and VT codes to obtain results where the spherical natures of the model clouds, 
e.g., limb-brightening and the geometrical size differences (see Section |4^ ), are not important by setting the 
virtual telescope beam size {0d) much smaller than the cloud size {2Rcioud)- This permits us, in effect, to 
obtain the line intensities from a plane-parallel cloud with a finite thickness whose front and back surfaces 
are exposed to external far-UV radiation. The surfaces of the cloud are perpendicular to the line-of-sight. 
These models are equivalent to setting 77 = 1 and = 1 in Equation ^: h = fi- The resulting [C 11] intensity 
is from both the front and back surfaces, because the [C 11] line emission is nearly optically thin (Stacey et 
al. 1991). The CO J = 1 ^ line is optically thick, so the resulting CO emission comes predominantly 
from the front surface. 

The VT code does not calculate the ro-vibrational lines of H2. We derive the H2 (1,0) S(l) emission 
from the column density, 7V(H2), which results from the plane-parallel PDR model. The electronically 
excited H2 molecules relax radiatively to the ground electronic state and then cascade through the 
vibrational energy levels by emitting ro-vibrational lines. In this process, the II2 line ratios are insensitive 
to nu and luvi and we can use a constant conversion factor between A^(H2) and II2 (1,0) S(l) intensity 
(Black & van Dishoeck 1987): 

jfront ^ 2.67 X 10"2^iV(H^) ergss-^cm-^sr-i , (24) 
where fl[°"^ is the intensity from the front surface of the cloud. 
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While the ro-vibrational hncs of H2 arc optically thin, the emission from the back surface, f^^^ , is 
affected by the extinction, Ak, through the cloud itself. Using Ak = 0.112Av (Rieke & Lebofsky 1985) and 
Equation we estimate the observed intensity from the back surface: 



log(/^T') = log(/i7"*) - ..Tm22" -2 • (25) 
^ 3.5 X lO^^cm 



As we will discuss in Section 5.4, the cloud size (in units of hydrogen column density N}j) has a lower limit 
to keep the CO molecules from being dissociated completely in an intense far-UV field {Ijjv > 10"^). Figure 
|l3| shows that pdustNu should be larger than 2 x 10^^ cm^^; therefore, 

fHT<^-^fHr'- (26) 

We will neglect the H2 emission from the back side of the cloud in the following discussion: fn^ — Z^™"*- 
In some cases, there may be a bit of additional flux from the back side which can penetrate through the 
thinner parts of the cloud. 

Figure O shows the [Cii], H2 (1,0) S(l), and CO J = 1 — > emission intensities from the PDR code 



and the MC/VT code for the two sided planar clouds. When dust absorption dominates (Equation 20), 
the H2 intensity increases as uh increases. On the other hand, when H2 self-shielding dominates, the H2 
intensity increases as Ijjv increases. In the LMC model (the right plot in Figure ^2|), the II2 intensity is 
enhanced by a factor of 10"-^ — lO^ "^ over the intensity in the Galactic model with the same nn and luv- 
The enhancement is mainly due to the different II2 self-shielding criteria in the Galactic model and in the 
LMC model. 



5.4. Model Emission Intensities Including the Effects of Spherical Geometry 

In order to understand the effect of limb-brightening and of differences in the physical sizes of the 
C^, II2, and CO zones in spherical clouds of varying metallicity, we ran the VT code, setting the virtual 
telescope beam size {Oed) to match the cloud size (2_Rc/o«d)- The resulting intensities are equivalent to 
setting r] = 1 in Equation |^: li ~ eifi. In Figure |l3|, we plot the model values of Icii, Ih2j and Ico versus 
cloud size. In Figures H, ||, and |lO|, we overlay the results on the observed data. Icii and Ico are obtained 
directly from the VT code, and the far-IR continuum emission and are from Equation ^ with cfir and 
eHn derived as follows: 



As we discussed in Section 4.3 , the effects of spherical geometry include both limb-brightening and size 
differences between the different emission regions in the cloud. The far-IR emission has very low optical 
depth and the II2 (1,0) S(l) emission is optically thin. The far-IR and II2 emission regions lie on the surface 
of the cloud and fill the telescope beam. We can analyze the spherical geometry effects for optically thin 
lines by projecting the three dimensional emission shell onto a two dimensional emission disk: 

AttR'^ dR 



Xi 2 / Xi 



where Xi is the depth from the cloud surface to the transition regions, e.g., Xq+ and Xh* which are defined 
in Section dT3; and Rdoud is the radius of the cloud. In the above equation, we assume that 71(112) /nn and 



7i(C+)/n/f are constant in the corresponding shells. For the far-IR emission, if we consider that the incident 
far-UV energy is conserved as the output far-IR energy (see Equation 0), we can write: 



^cloud 

We use Equation [2^ to obtain eH2^ ^'^^ Equation pi to obtain epij^^ = 2 



epiR = f^^^^ = 2 • (28) 



The CO line arises on the surface of the CO core. Because external UV radiation provides most of the 
heating, the gas temperature gradient is positive, dT/dR > (i? = at the cloud center). As we observe an 
optically thick line, the edge of the projected CO disk has a higher brightness temperature than the center, 
and we see limb-brightening. On the other hand, the size of the CO core (i?co) is smaller than the cloud 
size {Rcioud) which is affected by Ijjv and tih- When Rdoud is larger than X(j+^ and Ico is optically thick: 

V Rcioud ) V Rcioud J 

When we compare the spherical cloud models to plane-parallel models, the change in the CO core size 
(Rco) is more significant than the limb-brightening effect. We plot the CO intensity (by setting ?/ = 1) 
resulting from the numerical code in dotted lines in Figure This method also takes into account any 
additional effect of opacity variations on the local CO emissivity. As Rco ^ 0, the resulting Ico decreases 
rapidly. 



5.5. Applying the Model to the Data 



As we discussed in Sections |5.3| and |5.4|, the observed far-IR, [Cii], and II2 (1,0) S(l) line intensities 
depend almost exclusively on uh and luv, while the observed CO J = 1 ^ line intensity depends on nn, 
Ijjv, and Rdoud- In this section, we apply the model results to the data and estimate nn, Ijjv, and 
Rcioud in the regions we observed. 



5.5.1. Ih2 versus Icii 

In Table ||, the standard deviations of the observed \og{lH2/Icii) are smaller than those of other ratios 
(0.1 — 0.2 versus 0.2 — 0.4), because the values of fcii and ///^ depend similarly on nn and luv (see Figure 
p^ ), and in the case of spherical clouds, the [Cii] and II2 intensities are not very sensitive to variations in 
cloud size, i.e., 2/3 < (ec/ZjEffa) < 2 (see Equation |2^ and Figure p^) . 

In Figure ||, most of the data from NGC 2024 are within the model grids at 3.2 < lognn < 4.2, and 
1.5 < logljjv < 3. There are positions (at Aa — —17', —16', and —15') in NGC 2024, however, where the 
Ih2 intensities lie at least a factor of 2.5 above the model grids. At these positions, the clouds may be small 
enough to be transparent for II2 emission from the back side of the cloud (see Equation |2^) . Based on our 
experience with other positions, it is extremely unlikely that the II2 level populations at these positions are 
thermalized by effects present in high density PDR's (riH > 5 x 10'' cm^"^) or by shocks (see Section 
The ratios of \og{lH2 / Icii) in Orion are larger than those in NGC 2024. The Orion data in Figure ^ lie 
slightly below and slightly above models for logn^f ~ 4.7, implying that the gas density in Orion is higher 
than that in NGC 2024. 
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The mean ratios of log (//f 3 /Jc//) in 30 Doradus, N159, and N160 are similar to each other, suggesting 
that the observed clouds in the LMC have similar densities, and that coUisional de-excitation does not 
affect the H2 level populations (see Table ^ . In Figure ||, most of the data in the LMC match models in 
the range of 3.7 < logn/f < 4.7 and 2 < log/yy < 3. 

By comparing the observed data and the model results in Figure |[ we measured the median density 
and far-UV field for each source. With this derived nn and luv, we found the C'*' region depth, X(j+, using 



Figure 11, The values of tih, IuV: and X(j+ determined in this way are listed in Table 01 



5.5.2. Icii versus Ico 

Both cloud size differences and metallicity differences affect the observed Icii/Ico and the Ih^/Ico 
ratios. For the clouds in NGC 2024, the results of Jaffe et al. (1994) implied that the size of the clouds 



{Rcioud) in the western edge zone is smaller than in the cloud proper zone (see Section 11), while the C+ 



region depth, in the cloud is not much different. As we discussed in Section 5.4, for a given nu and 

luv, Ico is proportional to eco- In the LMC, the dust-to-gas ratio, p^^f , is lower than in the Galaxy by 



the factor of ^ 4, and the far-UV radiation penetrates deeper into the clouds (Figure |T^ ): 



vLMC ^GAL 
^ Pdust 
j^GAL — „LMC 



(30) 



C+ h'dust 

LMC ie locc tlior, ^GAL 



If we assume that the cloud sizes in the Galaxy and in the LMC are same, then e^.^ is less than ^co ■ 
^Co'~' — ^CO^' the cloud size in the LMC must be larger than the cloud size in the Galaxy (see Figure 14). 



The Icii /Icq ratio is also affected by the strength of the UV field. Mochizuki et al. (1994) compared 
the log(/c////co) ratios in the LMC with the ratios in the Galactic plane. Their log-ratio in the LMC, 
4.4 ± 0.3, agrees well with our ratio in 30 Doradus, N160, and N159 (4.2 — 4.8, see Table p[); however, their 
ratio in the Galactic plane, 3.1 ± 0.5, is much smaller than the ratio (3.7 — 4.3, see Table ||) in Orion and 
NGC 2024. We argue that the Galactic plane at —12° < I < +26° has a different physical environment from 
active star formation regions like Orion and NGC 2024. Most of the Iwiim/ Iwotim ratios at the positions in 
the Galactic plane are between 0.2 and 0.3 (see Figure 2 in Nakagawa et al. 1995). Applying these ratios to 
Equations |^, ^, and ^ (see also Equation ^) , we can estimate that the far-UV fields in the Galactic plane 
are in the range of 1.5 < \ogIuv < 1-7. The far-UV field is, therefore, much lower than toward the positions 
we observed in Orion and NGC 2024 (1.9 < logluv < 3.2). Our model results in Figure [T^ show that, at 
logljjv = 1-6 and logn^f — 3.7, the log(/c///^co) ratio is 3.2, which agrees with the observed value of the 
Galactic plane in Mochizuki et al. (1994). On the other hand, the far-UV fields in 30 Doradus, N160, and 
N159 are very high (2.0 < log/j/y < 3.8). Even the extended regions in the LMC in which Mochizuki et 
al. (1994) measured the log(/c;////co) ratio show higher far-UV fields (1.8 < log/f/y < 2.7, derived from 
the IeOfj.m/ IioOfim ratios) than those in the Galactic plane. The more intense far-UV fields, e.g., in Orion, 
NGC 2024, 30 Doradus, N160, and N159, dissociate more CO molecules and have deeper C"*" regions, X(j+, 
in Equation ^9]; therefore, the Icii/Ico ratio becomes larger. 

Figure^ plots the observed data, Icii/Ifir. versus Ico/Ifir, and overlays the results from the planar 
cloud models and the spherical cloud models. For the spherical cloud models, the tick marks along dashed 
lines mark the position in the space of models with different column densities, Nni'^Rcioud)- The locus of 
each model comes from the value of Icii—^ciifcii and Ico=^cofco determined for that column density. 
The solid lines and dotted lines show the results of the planar cloud model {Icii^fcii and Ico=fco)- 
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Only the data from the cloud proper zone in NGC 2024 are within the planar cloud model space, and most 
of the observed data deviate systematically toward lower Ico I Ifir or higher Icii /I fir from the model 



results. The models do not include spherical geometry effects (see Section x3), and Equation |27| implies 
that the observed Icii data should not be significantly affected by going from plane-parallel to spherical 
geometry. Therefore, the offset of observed points from the models on the Icii versus Ico plot in Figure || 
is solely a result of the effect of spherical geometry on the observed Ico emission (eco)- Given the observed 
value of Icii / Ifir, we measured how much the observed Ico /Ifir is shifted from the point expected from 
the models for a given density and far-UV intensity on the model grids in Figure |^: 

Alog(/co/-^_F/fl) = log eco • (31) 

Applying the measured Xc+ and eco to Equation |2^, we derived the column densities through the clouds 
{Nui'^Rcioud) in Table |l^) under the assumption that the clouds are spherical. For the data from the cloud 
proper zone in NGC 2024, we use the uncertainty in the observations to set the lower limit, i.e., eco ^ 

10-0-2. 

Madden et al. (1997) also estimated eco of the nearby irregular galaxy, IC 10, using far-IR, [Cii], and 
CO J = 1 — > data. In Table |l^, we list the adopted values of luv and eco at positions A and B in IC 
10. Using the dust-to-gas ratio of IC 10, p^^l^ — 0.2, and the FDR models, we estimate iV(C+) and then 

Nui'^Rcloud)- 



5.5.3. Ih2 versus Ico 

Models indicate that, like Icii, Ih2 is not affected strongly by changing from planar to spherical 
geometry. Therefore, the similar offset of the data from models in Figure |l^ is probably due to the 
inappropriateness of applying the plane-parallel models for predicting Ico- As were the data in Figure 
^, only the data from the cloud proper zone in NGC 2024 are within the model space, and most of the 
observed data deviate systematically leftward from the model results. 

The data from the SMC, which are not available in Figure ^, also deviate toward lower log{Ico / Ifir) 
or higher \og{I H2 1 1 f i r) from the LMC model results, if we assume that the SMC models without spherical 
geometry effects are not very different from the LMC models. In the SMC, the observed II2 emission could 
be higher than the model results as a result of the unreddened H2 emission from the back side of the cloud. 
Because [C 11] data are not available and the II2 data are not certain, we cannot measure e^^*^ from the Ih2 
versus Ico plane. We can still estimate a lower limit on the column density of clouds in the SMC, assuming 
eco > 10"^. Applying the dust temperature of 10^'^^ K from IRAS 60/im/lOO/im ratios (see Figure |^) to 
Equations ^ and ||, 

log /c/y - -5.65 + 5 log Tdust , (32) 

we obtain a far-UV field {luv) of lO^'® at the observed positions in the SMC. Assuming \ognH — 3.7, we 
get Nh{X^+"^) = 4.5 X 10^2 cm^^. Therefore, the central column density (Nh) of the clouds in the SMC 
should be at least 1.3 x 10^^ cm"2(see Table |l^). 



5.6. Cloud Column Densities in Galactic and Magellanic Complexes 



We estimate cloud column densities in the molecular complexes in the Galaxy and in the Magellanic 
Clouds. We first derive nn and luv by comparing the observed data and the model results in the IH2/IFIR 



- 20 - 



versus Icii/Ifir plot (Figure ^). We take a median value of the data in each object, except in N159 where 
we consider each position separately. Using the values of nn and Ijjv derived from Figurejq, we get Xq+ , 
the depth from the surface of the cloud to the C+-to-C-to-CO transition layer, in Figure We also get 
ccoi the parameter that accounts for spherical geometry, by comparing the data and models in Figure ^ 
When we do not have or H2 data, we assume that lognn = 3.7, and estimated Ijjv from dust color 
temperature using leofim/ IioOfj^m and Equations |^ and We finally derive the central column densities of 
clouds using Equation |2^, and list them in Table ^ (see also Figure [ij). Note that the column densities 
derived in this way are a measure of the number of molecules between the Hii/H boundary and the cloud 
center, not of the overall column denisyt through the GMC complexes. 

The cloud column density measured in the western edge zone in NGC 2024 is at least a factor of 5 
smaller than the column density measured in the cloud proper zone in NGC 2024. At the four data positions 
in N159, the variations of uh and Ijjv are not significant, e.g., 3.7 < lognn < 4.1, and 2.3 < logluv < 
2.7; however, the Ico intensities change by a factor of 10^'^, indicating that the column density varies more 
than a factor of 5 within the N159 region. The cloud sizes represented by the hydrogen column density for 
a given density, hh = 5 x 10^ cm'^, in the LMC {N^^'^^ = 3 - 15 x 10^^ cm'^), the SMC (> 1.3 x 10^3), 
and IC 10 (6 - 10 x 10^^ cm'^) are bigger than those in the Galaxy {N§^^ = 1 - 4 x 10^^ cm"^) 



5.7. Column Density and Metallicity 

Table ^ and Figures |l^ and ^ show that the column density and, for constant density, cloud size 
depend roughly linearly on metallicity and therefore the visual extinction through clouds is constant as one 
goes from the Galaxy to the LMC and SMC. 

McKee (1989) presented a theory of photoionization-regulated star formation (see also Bertoldi & 
McKee 1996). According to his theory, the rate of low-mass star formation is governed by ambipolar 
diffusion in magnetically supported clouds, and the ionization of gas is mostly by far-UV photons. He 
suggested that the mean extinction of molecular clouds in dynamical equilibrium is in the range of Ay = 
4 — 8 mag. Using Equation |l^, the range of extinction corresponds to: 

0.64 X 10^2 < pdust Nh < 1.3 x 10^2 cm-2 . (33) 

This range of column density roughly agrees with the empirical values of Galactic GMCs using Ico and 
virial theorem: Nh = 4 — 8 x 10^^ cm~^ (Scoville et al. 1987; Solomon et al. 1987). In regions with 
higher incident UV fields, the equilibrium range should move to higher column density. Our results from 
the Galaxy, the LMC, the SMC, and IC 10 are also within Ay =4—16 mag, and agree with the theory to 
within a factor of two for the range of observed incident fields. 



6. Conclusions 

We observed H2 (1,0) S(l), (2,1) S(l), and (5,3) 0(3) lines from the Magellanic Clouds, and detected 
emission even in regions where CO had not been seen during the ESO-SEST Key Program, implying that 
H2 is present in these regions. With deeper CO observations, we also detected weak (0.5 — 5 Kkms"^) CO 
emission in PDRs where the previous ESO-SEST Program did not detect CO J = 1 — > 0. We conclude 
that, in all regions where we could detect H2, CO is not completely dissociated and that the CO J = 1 ^ 
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emission is still optically thick in the Magellanic Clouds. Our observations may have consequences for 
studies using CO to trace molecular mass. 

Our spherical-shell cloud model demonstrates the effects of limb-brightening and differences in the 
physical sizes of the C"*", H2, and CO zones on the observed line strengths. The H2, [Cii], and CO line 
intensities from the plane-parallel cloud model are not very sensitive to the metallicity. When geometric 
effects are taken into account using spherical models, the H2, and [Cii] line intensities normalized to the 
far-lR intensity are also not very sensitive to the metallicity; however, the CO intensity is proportional 
to the surface area of the CO core. As we decrease the cloud size (or total column density, Nh), the CO 
intensity normalized to the far-IR intensity decreases. 

We compiled data of H2 (1,0) S(l), CO J = 1 — > 0, [Cii], and far-IR in the star formation regions 
in the Magellanic Clouds and in the Galaxy, and compared them with simulated line intensities from a 
plane-parallel cloud. The data in the I H2 1 1 fir versus Icii/Ifir plot agree with the model results, and 
show no significant difference between the Magellanic Clouds and the Galaxy. We use the Ir^/Ifir versus 
Icii/Ifir plot to estimate uh and luv at the positions we observed. In the Icii/Ifir versus Ico/Ifir 
plot and the IH2/IFIR, versus Ico/Ifir. plot, the data from the western edge zone in NGC 2024 and from 
the Magellanic Clouds are shifted to the lower Ico than the model results, which can be explained as 
an effect of spherical geometry. We measured the amount of the effect, and used the result to derive the 
column density through the clouds. 

On galaxy scales, the derived average column density and therefore size of the clouds appears to 
increase as the metallicity of the galaxy decreases. Our observational result, that the mean extinction of 
clouds is constant and independent of their metallicity, is consistent with the idea that photoionization may 
play a role in regulating star formation. 
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Table 1. H2 Lines and Instrument Parameters 



H2 Line 


A'' 


Date 


Bed'' 


ire 










/im 




arcsec 






kms" 


-1 


(5,3) 0(3) 


1.61308 


1995 Oet 


81.0 


23.0 


125 


104 


163 


(1,0) S(l) 


2.12125 


1995 Oct 


81.0 


25.4 


94 


125 


196 






1994 Dec 


88.9 


26.4 


94 


120 


189 


(2,1) S(l) 


2.24710 


1994 Dec 


88.9 


26.4 


89 


127 


199 



"Wavelength in air (Black & van Dishoeck 1987) 

''Beam size (equivalent disk). The beam profiles are in Figure^ 

'^Effective Finesse of the Fabry-Perot interferometer including the reflectivity, parallelism, surface quality, 
and incident angles 

Order of interference 

"Full-width at half maximum of the instrument profile of an extended source. 

^AVbw = J' I{y)d'V / Ip^ak- is the instrument profile of an extended source. 



Table 2. Frequency Switching Parameters 



Object 


H2 Line 




ON 




OFF 






lira 


km s^^ 


HTCi 


km 


SMC 


(5,3) 0(3) 


1.61386 


+124 


1.61506 


+347 




(1,0) S(l) 


2.12219 


+112 


2.12069 


-100 


LMC 


(5,3) 0(3) 


1.61464 


+274 


1.61584 


+477 




(1,0) S(l) 


2.12321 


+261 


2.12171 


+49 




(2,1) S(l) 


2.24935 


+249 


2.24785 


+49 
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Tabic' 3. Obj(x-t List 



Object'' Q:i950 <5i950 Reference 



DEM S 16 00 43 33.4 -73 39 05 1 

DEM S 37 00 46 16.2 - 73 32 47 1 

LI-SMC 36 00 44 50.5 -73 22 23 1 

30 Dor 05 39 11.5 -69 06 00 2 

N159/N160 05 40 18.2 -69 47 00 3 



(0, 0) Position 

References. — (1) Based on the CO survey of the ESO-SEST Key Program (Rubio et al. 
1993); (2) The CO peak (Poglitsch ct al. 1995); (3) Center of the CO map in the N159/N160 
region (Booth 1993; Johansson et al. 1994). Because the reference position of N160 is also 
N159 (0, 0) position, the center of N160 is close to the (-1, 7) position in this reference system. 
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Table 4. Intensities from the Magellanic Clouds 



Object 




A<5-'' 


H2 




CO 




CII 


fai-IR " 






Rcf. 












/ 


a 








CO;CII 




arcmin 


10"" 




K km s 


-1 


10-* 


10-" 


K 






DEM S 16 


+0.11 


-0.12 


3.43 


0.66 


2.19 


0.11 




0.43 


49 


2 






+0.11 


+0.88 


1.60 


0.66 


2.39 


0.11 




0.27 


42 


2 




DEM S 37 


+0.38 


-3.35 


0.99 


2.20 


















-0.02 


-0.05 


2.4 


1.1 


0.584 


0.073 




0.65 


45 


1 




LI-SMC 36 


-3 





1.33 


1.77 
























3.2 


1.3 


6.20 


0.060 




0.33 


42 


1 




30 Dor 





-1.5 


4.2 


2.7 


2.87 




2.6 


14 


45 


3 


3 




-1.5 





5.3 


2.2 


3.45 




4.2 


18 


51 


3 


3 










10.8 


0.93 


10.5 




8.1 


35 


75 


3 


3 




+ 1.5 





8.6 


2.4 


2.23 




3.7 


7.8 


44 


3 


3 







+ 1.5 


3.8 


2.9 


3.39 




2.9 


6.6 


42 


3 


3 




3 


-9 


0.94 


0.75 




3" 




1.4 


43 


4 









-6 


1.63 


0.78 


0.746 


0.052 




2.6 


43 


1 






-6 





-0.28 


0.53 




3" 




1.2 


56 


4 




N159 


-0.5 


— 5 


2.1 


1.5 


25.0 


0.14 




0.31 


34 


5 






+ 1 


— 5 


0.94 


0.90 


46.0 


0.14 




0.32 


35 


5 






-2.5 





0.5 


1.9 


29.4 


0.14 


0.77 


1.8 


31 


5 


6 




-1.5 





4.5 


1.4 


47.3 


0.14 


2.3 


7.1 


46 


5 


6 










3.23 


0.99 


7.70 


0.052 


2.5 


8.2 


42 


1 


6 




+1 


+ 1 


4.8 


1.3 


19.8 


0.052 


2.9 


8.7 


56 


1 


6 







+3 


1.80 


0.72 


0.634 


0.052 


1.0 


1.7 


46 


1 


6 


N160 





+4.5 


0.8 


1.2 




3" 


0.49 


1.2 


37 


5 


6 




-2.5 


+7 


1.95 


0.80 


1.59 


0.073 


0.77 


2.6 


36 


1 


6 




-1 


+7 


6.45 


0.91 


9.89 


0.073 


2.0 


15 


76 


1 


6 







+7 


2.46 


0.71 


5.82 


0.073 


1.7 


6.5 


54 


1 


6 




-2 


+8 


1.86 


0.64 


4.25 


0.073 


1.1 


2.9 


39 


1 


6 







+8 


2.18 


0.56 


6.76 


0.073 


1.1 


1.8 


34 


1 


6 




-3 


+9 


-0.21 


0.65 




3" 


0.68 


0.97 


35 


5 


6 




-1 


+9 


1.98 


0.59 


2.76 


0.073 


0.75 


1.1 


34 


1 


6 




-2 


+ 10 


1.72 


0.51 


1.64 


0.073 


0.55 


1.0 


46 


1 


6 







+ 10 


2.05 


0.52 


3.86 


0.073 


0.36 


0.75 


43 


1 


6 




-1 


+ 12 


-0.6 


1.3 


4.49 






0.71 


41 


5 





''Offset from the object (0, 0) position. 

''Measured H2 (1,0) S(l) intensity with the \u statistical uncertainty. 

'^Velocity integrated CO J = 1 ^ intensity, 1 = J Tmb d-v. 1 Kkms" -'^ = 1.57 X 10~^ ergss"-*^ cm"^ sr"^. 
'^[CII] 158 fim intensity. 

"IRAS Far-IR continuum intensity. See Equation ^ 

'Dust temperature deduced from the ratio, /eOfim/^ioOMm- See Equation |2| 
^In units of ergss"' cm~^ sr~'. 

''Upper Umit that hes beyond the lowest contour (Tmb = 3 Kkms"') in the published CO maps (Booth 1993; Johansson 
et al. 1994). See also Figures and K 



References. — (1) This work; (2)Rubio et al. 1993; (3) Poglitsch et al. 1995; (4) Booth et al. 1993; (5) Johansson et al 
1994; (6) Israel et al. 1996. 
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Table 5. H2 Line Ratios 



Object 


Aa 


AS 

arcmin 


Intensity'' 
10-s 


(2,1) S(l) 

Ratio'' 


Intensity* 

IQ-'' 


(5,3) 0(3) 

Ratio'' 


DEM S 16 


+0.11 


-0.12 






0.76 ± 0.34 


0.22 + 0.11 




+0.11 


+0.88 






0.27 + 0.36 


0.17 + 0.24 


30 Dor 








4.0 + 1.5 


0.37 + 0.14 


1.48 + 0.37 


0.138 + 0.036 




+1.5 









1.42 + 0.51 


0.166 + 0.075 


N159 


+1 


+1 






0.65 + 0.41 


0.135 + 0.091 


N160 


-1 


+7 


3.6 + 1.9 


0.56 + 0.30 


2.36 + 0.41 


0.366 + 0.082 



^Measured intensity with the la statistical uncertainty in units of ergss ^ cm ^ sr 
''Intensity ratio with respect to (1,0) S(l) line intensity. 



Table 6. Comparing with Other H2 (1,0) S(l) Data in N159 



Ref. 


Position'' 


Beam size'' 




Flux^ 


Intensity" 




arcmin 


arcsec 


10-8 sr 


10-14 


10-s 


1 


(+1.11, +0.85) 


6x6 


0.85 


1.7 


20 


2 


(+1.11, +0.85) 


(pis 


3.1 


5.2 


17 


3 


(+1.11, +0.85) 


10 X 21 


5.0 


8.5 


17 


4 


(+1.0, +1.0) 


4)81 


121 


58 


4.8 



Oflfset from the N159 (0, 0) position (see Table 3). (+1.11, +0.85) corresponds to "N159 Blob" (Heydari- 
Malayeri & Testor 1985). 

'' ^13 denotes Bed = 13". 

Solid angle of the beam, 
■l In units of ergss— l cm— 
" In units of ergss— 1 cm— sr— l. 

References. — (1) Israel & Koornneef 1992; (2) Krabbe et al. 1991; (3) Kawara, Nishida, Sz Taniguchi 
1988 (4) This work 
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Table 7. Intensities from the Galactic Clouds 



Object'' 


Aa 


AS 


Ha 




CO 


CII 


far-IR 




Ref.'' 










a" 


/ 








CO;CII 




arcmin 


10" 


6 


Kkms-i 


10-4 


10-2 


K 




Orion 


+8 





23 


6 


34.0 


9.9 


56 


55 


i;2 




+10 





28 


5 


29.0 


6.1 


20 


39 


i;2 




+12 





5.2 


1.6 


29.5 


4.3 


8.3 


36 


i;2 




+1.84 


+9.83 


27 


4 


60.7 


4.5 


41 


42 


1;3 


NGC 2024:Edge 


-17 





9 


1.9 


13.4 


5.6 


3.3 


36 


4;4 




-16 





12 


3 


18.6 


6.7 


4.7 


37 


4;4 




-15 





13 


2 


14.4 


8.6 


4.9 


37 


4;4 




-14 





3.3 


1.2 


6.70 


5.0 


4.3 


37 


4;4 




-13 





4.5 


1.8 


15.6 


4.7 


4.4 


35 


4;4 




-12 





6.4 


1.9 


51.0 


5.6 


6.6 


38 


4;4 




-11 





8 


2 


84.6 


6.3 


10 


41 


4;4 


NGC 2024:Prop 


-9 





9 


3 


110 


8.8 


19 


49 


4;4 




-8 





6 


3 


132 


11 


15 


43 


4;4 




-6 





8 


3 


154 


10 


26 


47 


4;4 




-5 





11 


4 


130 


9.9 


55 


52 


4;4 


Orion (0,0): o;i96o 


= 5''32' 


"^49?0, <5i95o = 


-05° 


25' 16" (6*1 


Ori C). 


NGC 2024 (0,0): 





5''39™14!0, <5i95o = -01°57'00". NGC 2024:Edge and NGC 2024:Prop have the same (0, 0) position. 
See the text for the differences between NGC 2024;Edge and NGC 2024:Prop. 

H2 data is from Luhman & Jaife 1996 

In units of ergss-^ cm-^ sr-^. 

References. — (1) Schloerb & Loren 1982; (2) Stacey et al. 1993; (3) Luhman et al. 1997; (4) Jaffe 
et al. 1994 

Note. — See the table footnotes in Table 4 for other columns. 
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Tal)k^ 8. Intc 


'iisity 


Ratios Ixrfwc 


Hm H2 


(1,0) S(l). [Cn 


|. and CO J 


= 1^0 




Object 


N". 




/Icii) 


Iog(7o// 


/Ico) 




/Ico) 






Avg.b 


S. D.' 


■= Avg. 


S. D. 


Avg. 


S. D. 


The Galaxy 


15 


-1.85 


0.28 


4.04 


0.37 


2.19 


0.47 


...Orion 


4 


-1.53 


0.31 


4.01 


0.25 


2.48 


0.32 


...NGC 2024:Prop 


4 


-2.08 


0.13 


3.68 


0.03 


1.61 


0.14 


...NGC 2024;Edge 


7 


-1.91 


0.15 


4.26 


0.37 


2.35 


0.38 


The LMC 


17 


-1.70 


0.17 


4.37 


0.44 


2.71 


0.40 


...30 Dor 


5 


-1.81 


0.11 


4.82 


0.14 


3.01 


0.23 


...N160 


8 


-1.59 


0.18 


4.18 


0.21 


2.59 


0.20 


...N159 


4 


-1.78 


0.07 


4.19 


0.65 


2.41 


0.62 


The SMC 


4 










2.89 


0.40 



Number of observed positions. 

^ Average of \og{Icii / Ih2)- lH2t Icii, and Ico denote hne intensities of H2 (1,0) S(l), [Cll] 158 /im, 
and CO J = 1 0. 

standard deviation (a measure of how widely values are dispersed from the average value) of 

^Og(Icil/lH2)- 



Table 9. Input parameters for the PDR Code 



Object 


10^2 cm-2 




Zo'' 


Pdust 


10-"s-i 


The Galaxy 


1 


1 


1 


1 


5 


The LMC 


4 


0.25 


0.25 


0.25 


5 


The SMC" 


10 


0.1 


0.33 


0.1 


5 



» H column density of the model cloud. Nh = N(}i) + 2iV(H2). 

^ Carbon and Oxygen abundances normalized to the Galactic values: [C]'''^^ = 1.6 X 10^^, and [O]''^^ 
= 5 X 10"**. Because the abundances are uncertain, we use simplified values in the model calculation. 

Dust-to-gas ratio normahzed to the Galactic value. [Av/A''if]''^^ = 6.29 x 10~^^ mag/cm^^ (Bohlin 
et al. 1983; Black & van Dishoeck 1987). 

An atomic hydrogen cosmic-ray ionizing frequency. 

" Only for a standard model with nn = 5 x 10^ cm~'^and luv = 10^- 



Table 10. Derived Cloud Column Densities 



Object 


Pdust^ 




lognH 


log luv 


N{Xc+r 




NHi^Rcloud)" 
1022cm-2 


The Galaxy 


1 


> 15 












...Galactic Plane' 






(3.7) 


1.6 


0.29 


> 0.63 


> 1.4 


...Orion 




4 


4.7 


2.5 


0.27 


0.44 


1.6 


...NGC 2024:Prop 




4 


3.7 


2.5 


0.45 


> 0.63 


> 4.3 


...NGC 2024;Edge 




7 


3.7 


1.8 


0.32 


0.083 


0.89 


The LMC 


0.25 


17 












...30 Dor 




5 


4.0 


2.7 


1.5 


0.048 


3.9 


...N160 




8 


4.2 


2.7 


1.5 


0.17 


5.0 


...N159(-1.5,0) 




1 


4.1 


2.7 


1.5 


> 0.63 


> 15 


...N159(+l,+l) 




1 


3.9 


2.7 


1.7 


0.36 


8.4 


...N159(0, 0) 




1 


3.7 


2.7 


1.8 


0.15 


5.8 


...N159(0,+3) 




1 


4.1 


2.3 


1.2 


0.063 


3.2 


The SMG' 


0.1 


4 


(3.7) 


2.6 


4.5 


> 0.1 


> 13 


IC lO'^'S 


0.2 


2 












...Position A 




1 


(3.7) 


2.7 


2.2 


0.3 


9.9 


...Position B 




1 


(3.7) 


2.6 


2.0 


0.1 


5.8 



Dust-to-gas ratio normalized to the Galactic value. 
^ Number of observed positions. 

" Depth from the surface of the cloud to the C+-to-C-CO transition layers. 

Parameter of the spherical geometry effect. 

" The central column density of the spherical cloud. Note that Nni'^Rcloud) = "^Rdoudi^H, and Rdoud 
X(-,+ + Rco- 

' We assume that lognn = 3.7. We derive log Ijjv from IRAS far-IR data. See the text for details, 
s We adopt results from Madden et al. (1997). 
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-60 60 -60 60 

offset [arcsec] 

Fig. 1. — Beam profiles of the University of Texas Near-Infrared Fabry-Perot Spectrometer at the CTIO 
1.5 m telescope in 1994 December and 1995 October. 
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Fig. 2. — Spectra of H2 (1,0) S(l) and (5,3) 0(3) emission lines (solid histogram) at the (0,0) position of 30 
Doradus in units of 10~^ ergss"^ cm~^ sr~^ (kms~^)~^. The smooth curve shows a fit to the spectra with 
the instrument parameters using GaussFit (Jefferys 1990). The overlapped telluric OH lines (dotted line), 
measured at an off-source position, are scaled down by 20 on the H2 (1,0) S(l) spectrum and by 100 on the 
H2 (5,3) 0(3) spectrum. 
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4 2 -2 -4 -6 4 2 -2 -4 -6 

Aa['] (5^ 39"" 11'. 5) 

Fig. 3. — Intensity maps of the 30 Doradus region. The (0, 0) position is 01950 — 5'^ 39'" 11!5 and (S1950 
= —69° 6' 0". The axes show the R.A. and Dec. offset in arcmin. Left: IRAS far-IR continuum map 
(thick hnes) and new CO J = 1 ^ map (thin lines near (0, -6')). The far-IR intensities are calculated from 
the observed leofim and Iioo^m using Equation |l|. The CO intensity data are obtained from this work, and 
the crosses show the observed positions. The contours are spaced at logarithmic intervals: Ifir = 0.398, 
0.631, 1.00, 1.58, 2.51, 3.98, 6.31, 10.0, 15.8, 25.1, 39.8 (i.e., 10°-2", n =-2, -1, 8) in units of 10^^ 
ergss"icm-2sr-i; and Tmb = 0.398, 0.631, 1.00, 1.58 (i.e., lO'' ^", n = -2, -1, 0, 1) in units of Kkms^i. 
Right: Plotted H2 (1,0) S(l) data (within the circles) overlaid onto the published CO J = 1 ^ map. The 
numbers are in units of 10~^ ergs s~^ cm~^ sr~^, and the size of the circles are the equivalent-disk size of the 
UT FPS beam. For the CO data (Booth 1993, Johansson 1994) the contours are spaced linearly: Tmb = 3, 
6, 9, ... in units of Kknis"^. 
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2 -2 -4 2 -2 -4 2 -2 -4 

Aa['] (5^ 40"" 18.2') 

Fig. 4. — Intensity maps of the N159/N160 region. Left: Far-IR continuum map. The contour intervals 
are spaced at logarithmic intervals: IpiR = 0.158, 0.251, 0.398, 0.631, 1.00, 1.58, 2.51, 3.98, 6.31, 10.0, 15.8 
in units of 10"^ ergss^"'^ cm^^ sr^^. The far-IR peak at (— 1', +7') is near N160 and the peak at (— 1',0) is 
near N159. Middle: New CO J = 1 ^ map. See the caption in Figure || which shows the same data. 
Right: Plotted H2 (1,0) S(l) data (within the circles) overlaid onto the published CO J = 1 ^ map. See 
the caption for the right map of Figure H. 
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CO in N159/N160 



1 — ' — ' — ' — \ — ' — ' — ' — \ — ' — ' — r 




I \ ^ I > I > ^ I \ ^ 1 ^ I 

2 0-2-4 

Aa['] (5^ 40™ 18.2') 

Fig. 5. — New CO J = 1 ^ intensity map of the N159/N160 region. The contours are spaced at logarithmic 
intervals: (for the (-l',+8.'5) region) Tmb = 0.398, 0.631, 1.00, 1.58, 2.51, 3.98, 6.31, 10.0, 15.8 Kkms"!; 
(for the (0,+3') region) 0.398, 0.631, 1.00, 1.58 Kkms'^ (for the (+0.'5,+0f5) region) 0.398, 0.631, 1.00, 
1.58, 2.51, 3.98, 6.31, 10.0, 15.8, 25.1 Kkms~^. The observed positions are plotted in plus signs. 
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Fig. 6. — [Cii] 158 /iin intensity versus CO J=1^0 intensity in units of ergss^^ cm^^ sr^^. The left plot 
is replicated from Figure 4 in Jaffe et al. (1994). The right plot shows the distribution of the Galactic 
data which we are using in this work. We also overlay the results from the plane-parallel models, using the 
Galactic parameters. See Figures ^ and || for explanations. 
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Fig. 7. — Far-IR intensities (Equation |^ and Table-94 an6 7) versus Tdust (Equation || and Tables 4 and 
7). IpiR is in units of ergss"^ cni~^ sr~^. The data in the left plot are from the Galactic clouds and, in the 
right plot, from the LMC. The dotted line shows the model in Equation |[ We do not include the effects of 
spherical geometry {epiR, see Section 5.4) in the model calculation. 
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Fig. 8. — I H2 1 1 FIR versus Icii / Ifir- The data in the left plot is from the Galactic clouds and the right 
plot from the LMC. We overlay the results from the two-sided plane-parallel models (see Section |5^). We 
used the Galactic parameters for the model results in the left plot and the LMC parameters in the right plot 
(see Table 9 for the list of parameters). In the model grids, the solid- lines and the dotted lines connect the 
same nn and the same luv respectively. n3.7 denotes nn — 5 x lO'^cm"'^, and i3 denotes Ijjv = 10^. 



- 39 - 




Fig. 9. — Icii/Ifir. versus Ico/Ifir- See Figure ^ for the explanations of the model grids of the solid 
lines and the dotted lines. We also plot models which include the effects of spherical geometry (shown as 
dashed lines, see Section 5.4). The dashed lines connect the spherical models with the same nn and luv 
with different cloud sizes: the upper dashed line is for uh = 5 x 10"^ cm^'^ and Ijjv = 10^; and the lower 
line for nn = 5 x 10"^ cm"-^ and luv — 10'^. The plus sign on the dashed line marks the cloud size spaced 
by 0.1 in logarithmic-scale. iV22.3 denotes Nni^Rcioud) = '^nnRcioud = 10^^ '^ cm^^. Also see Figure for 
another plot of these data. 
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1 2 3 4 1 2 3 4 



log (luv) 

Fig. 11. — luv versus the depths {Nh{X) = Xuh) from the surface of the cloud to the H-to-H2 transition 
layer (where the abundance of H becomes same as that of H2, shown as solid lines), and to the C^-to-C-to- 
CO transition layers (where the abundance of C+ becomes same as that of CO, shown as dashed lines). We 
changed the initial parameters, e.g., nn, luv , and Z, for each model. n2.7 denotes nn = 5 x 10^ cm~^. At 
fixed density, the dashed lines and the solid lines show the change of X(j+ and Xh* respectively. The dotted 
line divides the {Iuv^ t^h, Xh*) space into the H2 self-shielding dominant space and the dust absorption 
dominant space (see Equation |2^) . 
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Fig. 12. — luv versus [C ii] (dashed line), H2 (1,0) S(l) (solid line), and CO J = 1 — > (dotted line) emission 
line intensities at fixed density. We assume that 77 = 1 and the intensities do not include the geometry effect 
(ci = 1): Ii = MnH,Iuv,Z). 
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Fig. 13.— [Cii] (dashed line), H2 (1,0) S(l) (solid line), and CO J = 1 -> (dotted line) line intensities 
versus Rdoud at {uh = 5 x 10^ cm~^; Iuv= 10^, 10'^). 12 and j3 denote luv = 10^ and 10'^ respectively. 
We assume that rj — 1 and the intensities include the effect of spherical geometry: li — eifi{nH, luv , Z). 
Icii and Ico are directly from the VT code, and is from e//2/i?2 with Equations 24 and E^. We plot 



only the range of Nh where significant changes take place in each model. The upper axis shows the cloud 
radius in pc, and the lower axis the central column density of the spherical cloud in cm^^: Nni^Rcioud) = 
'^Rcioud.nH- Note that Nh(X) in Figure |ll| is Xuh, where Xc+ = Rdoud - Rco- 
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The LMC 




Fig. 14. — Schematic of chemical structures and sizes of typical clouds in the Galaxy and the LMC. Assuming 
that the far-UV field (luv) and the gas density (nn) are same in the clouds in both galaxies, the cloud size 
in the LMC is four times bigger than that in the Galaxy, i.e., the typical size of the star-forming clouds is 
inversely proportional to the metallicity of the galaxy. 
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Fig. 15. — Central column density of spherical cloud versus dust-to-gas ratio from the Galaxy, the LMC, 
the SMC, and IC 10 (see also Table ^0|). Note that the cloud size is presented by the hydrogen column 
density for a given density, nn = 5 x 10"^ cm~'^: Nni'^Rcioud) — '^Rcioudnn, and pdust is normalized to the 
Galactic value. The area of the filled circle represents the number of observed positions. The arrow shows 
the lower limit. To compare with the theory of photoionization-regulated star formation, we overlay the 
constant optical depths at Ay = 4, 8, and 16 mag. 



